Metamagnetism and Lifshitz Transitions in Models for Heavy Fermions. 



M. Bercx and F. F. Assaad 
Institut fur Theoretische Physik und Astrophysik, 
Universitat Wiirzburg, Am Hubland, D-97074 Wurzburg, Germany 

We investigate metamagnetic transitions in models for heavy fermions by considering the doped 
Kondo lattice model in two dimensions. Results are obtained within the framework of dynamical 
mean field and dynamical cluster approximations. Universal magnetization curves for different 
temperatures and Kondo couplings develop upon scaling with the lattice coherence temperature. 
Furthermore, the coupling of the local moments to the magnetic field is varied to take into account 
the different Lande factors of localized and itinerant electrons. The competition between the lattice 
coherence scale and the Zeeman energy scale allows for two interpretations of the metamagnetism 
in heavy fermions: Kondo breakdown or Lifshitz transitions. By tracking the single-particle residue 
through the transition, we can uniquely conclude in favor of the Lifshitz transition scenario. In this 
scenario, a quasiparticle band drops below the Fermi energy which leads to a change in topology of 
the Fermi surface. 
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I. INTRODUCTION 

Kondo lattice systems are states of matter whose low 
temperature macroscopic properties are dominated by 
strong correlations between Bloch fermion states and lo- 
cal spin moments. They can host various, sometimes 
competing orders and are therefore susceptible to tun- 
ing of external parameters. How the strongly entangled 
Kondo state evolves when competing mechanisms ap- 
pear, constitutes a vibrant area of research!^] Prototyp- 
ical heavy fermion materials are intermetallic compounds 
with the rare-earth elements Ce or Yb that deliver almost 
localized /-electrons. 

When an external magnetic field is applied to cer- 
tain fermion systems, unexpected non-linear behavior of 
the mag netiz ation at a well-defined field value enters 
the stageP^l Equally, distinct anomalies of thermody- 
namic quantities and in tran sport measurements occur 
at the same magnetic fieldP^ This phenomenon has been 
dubbed metamagnetism. At the critical fi eld, the heavy 
electron Fermi surface changes its topologyPSfini Recent 
experiments witness a pronounced, first order metamag- 
netic transition (MMT) in the heavy-fermion paramag- 
net CeTiGe. 5 Metamagnetism has been known to occur 
m CeRu 2 S# and, amongst other fermionic systems, a 
pressure-tuned first order MMT has been observed in bi- 
layer ruthenatesP^ 

The thermodynamic signatures of heavy fermion com- 
pounds have been related to a metamagnetic quan tum 
critical endpoint of the Ising universality class J12U2I 
MMT in heavy fermion syste ms h as been addressed by 
static mean- field (MF) studies^^ 5 - 1 - presupposing a con- 
tinuous transition - and by dynamical mean field theory 
(DMFT). The magnetization profile in Kondo systems 
has been shown to be closely related to the quasiparti- 
cle coherence.^ Also, crystal field effects have been in- 
cluded in a DMFT studyP^l Apart from heavy fermion 
systems, the metamagnetism of itinerant electrons has 
been addressed by MF methods^, functional renormal- 



ization groujPl and DMFTP 

Lifshitz transitions are quantum phase transitions which 
invoke a topological change of the Fermi surf acel ELH221 Lif- 
shitz transitions and Kondo breakdown scenarios have 
been investigated in fermionic large-TV approaches! 24 * 25 ! 
Zeeman-driven Lifshitz transitions were shown to explain 
many anomalies in thermodynamic and transport mea- 
surements of certain heavy fermion metals! 2 ^ 
This study is motivated by the interplay of two com- 
peting energy scales, the lattice coherence scale and a 
magnetic Zeeman scale. By varying the magnitude of the 
Lande factors we can show that the metamagnetic transi- 
tion occurs when both scales are comparable, thus allow- 
ing for interpretations based on Kondo breakdown or Lif- 
shitz transitions. The single-particle residue is measured 
as a function of magnetic field throughout the MMT and 
is shown to be consistent with the picture of a coherent 
band dropping below the Fermi energy at the transition. 
We supplement our analysis by single-particle spectral 
data. Our results clearly point towards Lifshitz physics as 
the key player in the MMT in models of heavy fermions. 
We draw this conclusion based on a dynamical Cluster 
Approximation (DCA) calculation of the Kondo lattice 
model with a Hirsch-Fye quantum Monte Carlo solver. 
The paper is organized as follows. Section [IT] introduces 
the model Hamiltonian and Sec. lIIII reviews the DCA im- 
plementation. Sections |IV| and M contain the results of 



this study. We finish with a discussion (Sec. VI I and the 
conclusion (Sec. VII I. 



II. MODEL 

The essential aspects of heavy fermion sy stems are 
captured by the Kondo lattice model (KLM)P323 The 
KLM is an effective low-energy model which is obtained 
upon integrating out the valence fluctuations of the /- 
orbitals in the periodic Anderson modelP^Hl j n particu- 
lar, the model captures the crossover from independent 
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magnetic impurities embedded in a metallic host to a co- 
herent heavy fermion state. The KLM at half-filling has a 
unique spin singlet, insulating ground state^that is adi- 
abatically connected to the trivial band insulator of the 
non-interacting periodic Anderson modelP^ The weakly 
doped KLM exhibits a Fermi liquid ground stated 
We invest igate this model by means of DMFT and 
with a quantum Monte Carlo cluster solver. 
The cluster approximation is on spatial correlations 
which are essentially cut off by the cluster dimension. 
Temporal correlations that drive the Kondo effect are 
fully accounted for in DMFT and its cluster extensions. 
We study the KLM supplemented with Zeeman terms on 
the two dimensional square lattice, rl = rl t + rlj + ~Hb'- 



M = -i E (c«W + H.c.) + JX) S i- S i 

<i,j>,a i 

-fi B Bj2i9cS c Zii + 9 fS f Zti ). (1) 



The magnetic moments of itinerant (c) and local (/) or- 
bitals along the direction of the applied field are given by 
= Hb9cjS^,'1 . The couplings g c and gj are under- 
stood as parameters. Physically, this is motivated by the 
pseudo-spin nature of S[ f the spin degree of freedom 
originates from a Kramer doublet and can take large val- 
ues which in turn renormalizes the g- factor g/P^ 
Drawing on the rece ntly obt ained phase diagram of the 
two dimensional KL M 32 ' 35 ' 36 -! we concentrate on the para- 
magnetic side of the transition and consider a metallic 
state with the conduction band filling n c = 0.9. 
The lattice of impurities introduces the coherence scale 
T co h as the natural energy scaled The single-impurity 
Kondo scale, itself being the natural scale in a single- 
impurity model, is a local scaled 

Guidelines for these two scales in the KLM are pro- 
vided by large- N calculations.^ In the weak coupling 
limit (J/W -C 1, W is the bandwidth) and at small 
deviation from half- filling (1 — n c <C 1) a scaling of 
T co h oc T K oc We~ Poi - ei?)/J is obtained (poOf) is the free 
density of states at the Fermi level). 
Two dimensional Kondo systems are realized in su rface 
alloys, e.g. in the heavy fermion compound CePtsPSffl 
In the case that the distance to a continuous quantum 
critical point is large enough so that the zero dimensional 
Kondo effect dominates over spatial fluctuations one can 
expect that a similar scenario of competing energy scales 
applies to the three dimensional case. 
In the model (Eq. [I]) spin-orbit coupling is neglected 
which would generally cause the g-factor to be a tensor. 
Realistic modeling of heavy fermion materials requires 
a more sophisticated approach capturing these material 
specific features. Instead, the used model serves the pur- 
pose of describing the generic interplay between the mag- 
netic Zeeman scale and the coherence scale of the Kondo 
system which can lead either to the Kondo breakdown or 
the Lifshitz transition scenario. 



III. METHOD 

We use the Hirsch-Fye QMC technique to solve the 
KLM on small clusters, that contain two orbitals (DMFT 
limit) and four orbitals, respectively. Cluster approxima- 
tion schemes are particularly well designed to capture the 
Kondo physics since temporal correlations can be treated 
exactly The approximation is on spatial correlations 
that are short ranged in the present situation. The DCA 
is a fully causal, non-perturbative method which is sys- 
tematically improved by increasing the cluster sizes 1221311 
In the following, we outline our implementation for the 
KLM. 

A static mean-field perspective can provide insight into 
the low energy properties of the KLM.E3. It roots on the 
saddle-point approximation which is the exact solution 
of the SU(N) KLM in the limit of N = oo. However, 
it exhibits an unphysical phase transition instead of the 
Kondo crossover. Appropriate choices of the magnetic 
matrix elements that couple the impurity /-orbitals to 
the external magnetic filed can recover the smooth Kondo 
crossover even in the large- N limit of the KLM with an 
external magnetic field term.^ 

In order to solve the KLM, we implement the following 
HamiltoniarPE2l 



H -7/uH-?/, -jYl 



5>V 



fi.a + fi, a a 



(2) 



Here, W = % t + Hb and the Hubbard term Hu = 
if Ylii [ Ylia n l(T ~ l] " ^ as been introduced. Local spin- 
operators S{ are as usually mapped to auxiliary lattice 

fermions, S{ = \ Y. a ,f} f} a a oipfi,p- Their single occu- 
pancy is guaranteed for Uf — > oo and in this limit, the 
Hamiltonian |2]) is equivalent to the KLM (JlJ. 
The discretization j3 — MAt on the interval [0, /3] gives 
the partition function Z = Za t + C[(^ r ) 2 ]j with 



M f 

Zat =TrJ|| exp [ — At Ho] 

i=i 



D[A] exp 



(3) 

iArJXJX-*)" 

i a J 

x |d[0] exp - ArjJ2 (fu ~ ^z2( c L^ + H ' c -)) } 
= h[\,(f>] exp[-S eS [X,4>}] ■ 

In Eq. |3j the two successive Hubbard-Stratonovich (HS) 
transformation reduce the quartic fermion terms to 
quadratic terms. The integration measures D[A,</>] de- 
note integration over spatial and time indices of the fields 
and contain normalization factors. 

The saddle-point of the above defined action fulfills 
dSes/d<fin = dS c e/dXu = 0. Static mean-field theory is 
obtained by dropping the r-dependence in the HS fields, 
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and one can furthermore request the homogeneous solu- 
tion: (fin = 4>q ,Xu = Xq. The saddle-point equations, 



1 = ^S"Umf , 



(4) 

are then solved self-consistently. The respective mean 
field results for total magnetization and quasiparticle 
residues are discussed in Sec. llVl 

In order to go beyond mean field, a systematic 1/N 
expansion around the mean field solution can be 
performed.^ Instead, we integrate over all the field con- 
figurations through application of the HF-QMC algo- 
rithm. The trace in Eq. [3] can be carried out and ex- 
pressed as a determinant of the Green's function matrix 
g a . Then, the partition function 



Z 



At 



DfA, 



n det t^ 



(5) 



is sampled stochastically. In the actual implementation, 
two discrete HS transformation are usedPSl The Green's 
function matrix is measured and updated according to 
the Hirsch-Fye algorithm.^ During the simulation of 
Eq. [2j double occupancy of the /-orbitals can be sup- 
pressed to the desired accuracy. We take Ar = 0.25 
during the simulations. We have checked that smaller 
values of Ar do not alter the results. 
The cluster approximation amounts to considering the 
interaction terms of the Hamiltonian only on a subset 
M. of the lattice with N c sites, which naturally defines 
the extent to which spatial correlations are captured. We 
therefore solve the model 



H — iio + J ^ s c R ■ S R 

REM 



(6) 



by using the auxiliary Hamiltonian ([2j). Hq denotes the 
bath which is determined self-consistently. 
The DCA is naturally described in momentum space 
since it relies on coarse-graining of momentum space. 
Since the interaction part of the Hamiltonian is local, it 
is not affected by the coarse-graining. The model Hamil- 
tonian is solved on a finite cluster of N c sites that is 
embedded in a bath of N sites (N ^> N c ). Since N is not 
a limiting factor one can work directly in the thermody- 
namic limit. Therefore, the DCA interpolates between 
two limiting cases: the DMFT (N c — 1) and the finite 
lattice (N c = N). 

The DCA lattice self energy is a step-function in recip- 
rocal space: 

EE a C t f(K, w ) = ^^S Latt (K + k, w ), (7) 

k 

^Um v SE a C t f(K,a;) = E Latt (K ) c t ;). 

The step size is AK = 2ir/N c , the cluster momenta K 
define the centers of N c reciprocal cells and k denotes 



the k-points that lie within these cells. The DCA self- 
consistent scheme operates on the single-particle level of 
the self energies and it demands that E§j^ t01 .(K, oj) — 
£l^(K,u;). The self-consistent equations are: 



DCA 
Cluster 



(K,w) 



(8) 
-l 



= G^tvXK^)- 1 + EE a C t f (K, W ) - Ggpi tor (K, w ) 

Here, an effective bare Green function has been defined 

as 

g^ZLtK,^- 1 = GEStvXK,^- 1 + x™£ (k, w ) . 

(9) 

The DCA lattice averaged Green functions are coarse- 
grained averages over cell momenta: 



✓-.DCA 
^Latt^av 



(K,w) 



1 



JV^ w _ e (K + k) + / i-E£SA(K 



(10) 

The Green function Cluster Q^i w ) * s tne bare Green 
function that is the input for the cluster calculation. 
The cluster calculation yields the cluster Green functions 
G§^ cr (K,o;) that enter Eq. [8] Once the self-energy is 
converged, the DCA lattice Green function is computed: 

Gg£*(k, W )- 1 =c- e (k) + M -£Ea C u (M(k), W ) . (11) 

The function M : k — > K uniquely maps momenta to the 
reciprocal cells. 

The required CPU time of the HF-QMC algorithm scales 
as (/3A C ) 3 . 



IV. RESULTS 

The magnetic field tunes the interacting Kondo sys- 
tem (Eq. [I]) from strong coupling at low fields to weak 
coupling at high fields, i?/T co h 3> 1. This limit is adi- 
abatically connected to two copies of non-interacting c- 
electrons, spin split by the Zeeman energy, and fully po- 
larized /-moments. At low values of the magnetic field, 
the hybridized band is expected to shift in a rigid man- 
ner. At an intermediate energy scale, B ~ Tech, two 
different scenarios are conceivable: (1), a breakdown of 
the Kondo effect itself at the relevant energy scale or, 
(2), a continuous transition that preserves the quasipar- 
ticles. In scenario (1), the quasiparticle itself is destroyed 
by the magnetic field. The single-particle residue quan- 
tifies the overlap of the interacting wave function with a 
bare conduction electron wave function. Therefore, the 
loss of quasiparticle coherence has to manifest itself as 
a sudden drop in this quantity for both both spin pro- 
jections. In scenario (2), quasiparticles remain intact at 
the Fermi level. The spin dependent Fermi surfaces un- 
dergo Lifshitz transitions which modify their topology. 
As shown below, data for the single-particle residue and 
single-particle spectral function across the MMT support 
scenario (2). 

The mean-field solution, derived from the saddle-point 
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of the SU(N) KLM, can be seen as the best approxima- 
tion in quadratic fermionic terms to the fully correlated 
model. Therefore, in the case that DMFT/DCA calcu- 
lations support the notion of quasiparticles, the MF per- 
spective is legitimate. The MF results are intended to 
complete the above described scenario of Lifshitz transi- 
tions. 

The magnetization profile of a heavy fermion model 
system can serve directly as a measure of coher- 
ence. The plateau of the occupation number difference, 
m = o(n% + directly relates to the hybridization 
gap in the quasiparticle bandsP^ Its position is fixed 
to x = 1 — n c by the Luttinger sum rule. The physi- 
cal magnetization M = — |g = J2 a a(g c n^ + gfnfy does 
not generally display a plateau when the orbital couplings 
are not the same. 



A. Data collapse M (g f /g c , T/T coh , B/T coh ) 

For temperatures below the coherence scale T co h quasi- 
particle bands are formed via coherent superposition of 
the screening clouds of the local spins. To verify that the 
coherence scale is the unique underlying scale a data col- 
lapse of DMFT and DCA data is carried out by scaling 
the magnetization with T co h- Because the plateau width 
of the occupation number is a measure of the hybridiza- 
tion gap, a good estimate of T coh is obtained from the 
position of the second kink at B = i?L2, determined by 
the intersections of linear fits at fit = 100 and gf/g c = 1- 
The scaling then becomes: 

M {g f /g c , T, B, J) -> M {g f /g c , T/T coh , B/T coh ) . (12) 

Effectively, the Kondo coupling J/t has disappeared as a 
parameter in M. The data collapse is evident, as shown 
in Fig. [T] 

For all values of gf/g c the magnetization shows two 
pronounced kinks at B = -Bli.2- The driving mechanism 
that shapes the magnetization is rooted in the competi- 
tion of two energy scales: the dominant magnetic en- 
ergy scale gjusB and the Kondo scale T co h. At the 
second kink both scales become comparable, such that 
-Sl2 <x gj ■ The position of the second kink in depen- 
dence of the coupling ratio is shown in Fig. [2] and the 
data are in good agreement with the above argument. 
Increased Zccman coupling to the local spins provokes 
the intermediate, plateau-like region to decrease and ren- 
ders the increase at B = £>L2 much steeper. 
Static MF calculations succeed in reproducing the qual- 
itative shape of M (Fig. |3|. In the MF picture, the two 
kinks in the magnetization correspond to two Lifshitz 
transitions. 

At this point, the data collapse of the magnetization can 
be compared to a scaling approach of the resistivity in a 
recent cluster DMFT (CDMFT) study of the Anderson 
lattice model close to half-filling of the conduction band, 
which equally reveals the lattice coherence temperature 
as the single underlying energy scaled 
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FIG. 1. (Color online) T co h-scaled magnetization 

M(g f / g c ,T/T coh , B /T coh ) for (a) two-orbital DMFT and (b) 
four-orbital DCA. 

The respective coherence temperatures are: (a) T co h(gf/g c = 
1, J/t = 1.6) = 0.094t, T coh {g f /g c = 1, J/t = 1.4) = 0.057£ 
and T coh (g f /g c = 1, J/t = 1.3) = 0.044f; (b) T coh {g f /g c = 
1, J/t = 1.6) = O.lOt and T cob (g f /g c = 1, J/t = 2.0) = 0.181 



Our calculated metamagnetic curves, as shown in Fig. [T] 
bear notable similarity with recent experimental data of 
the paramagnetic heavy-fermion system CeTiGeP This 
is discussed in Sec. I VI I 



B. Single-particle quantities: residue Z a and 
spectral function y4 CT (k,tj) 

The analysis of the single-particle quantities is based 
on the observation, that the KLM has a Fermi liquid 
ground state for the chosen value of conduction band fill- 
ing, n c = 0.9.^1, and for zero external magnetic field. The 
calculations were performed for J/t — 1.3, g//g c = 4, 
fJt = 100 and f3t = 200. For these parameters, we identify 
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FIG. 2. (Color online) The magnetic field value of the second 
Lifshitz transition, _B L2 (g//<7 c ), agrees well with _B L2 (<?//<7c) oc 



FIG. 4. (Color online) The imaginary part of Matsubara self 
energies ImS^^J (iui n ) at values of the magnetic field close 
to B = Bl2- 




N 



FIG. 3. (Color online) Static MF results for the magnetization 
M and single-particle residue Z a . 



two Lifshitz transitions that occur at /xb-BliA ~ 0.002 
and at ^Buijt « 0.01075. The latter corresponds to 
the MMT. 

A Fermi liquid signature is the analyticity of the re- 
tarded self energy S(k, u>) around the Fermi energy such 
that E(k, lo) allows for polynomial expansion. Then, 
the single-particle residue reads, expressed with the k- 
independent Matsubara self-energy T^ MFT (iui n ): 



7DMFT 



= lim 

T-s-0 



1 



ImS™ T K) 



(13) 



The quantity ImS^ 1 (oj n ) across the MMT at B = B h2 
is displayed in Fig. [4j Evidently, the imaginary part of 
the Matsubara self energy is free of divergences for both 
spin projections at low frequencies w n . We take this as 
evidence for the continuous transition scenario. 
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FIG. 5. (Color online) Single-particle spectral function 
A a (k,uj) across the MMT at the lowest temperature, fit — 200 
(magnetic field values are indicated by arrows in Fig. [6]). The 
narrow distribution of spectral weight close to the Fermi en- 
ergy (dashed line) indicates that Kondo coherence remains 
across the MMT. 



The excitations are tracked by the single-particle spectral 
function, 



1. 



A CT (k,w) = — ImG£ att (k,w) , 

7T 



(14) 



see Fig. [5j The analytic continuation from imaginary 
time dependent QMC data has been performed with the 
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FIG. 6. (Color online) Magnetization M and single-particle 
residue Z° MFT across the MMT. (The arrows refer to the 
single-particle spectra of Fig. [5]) 



stochastic maximum entropy method.^ 
The single-particle residues Z a at fit = 100 and f3t = 
200 across the MMT are shown in Fig. ^b). Zi es- 
sentially follows the magnetization M(B) (Fig. |6[a)). 
Ai(k,uj) displays well defined quasiparticle weight across 
the MMT (Fig. [5]) and hence accounts for a metallic state. 

vanishes for an intermediate magnetic field range, 
close to Bli < B < B^- In this locked phase, no 
up-spin Fermi surface is present. At B = B^2, a topo- 
logical change of the Fermi surface occurs since one up- 
spin band crosses the Fermi level at the gamma point, 
(k x , k y ) = (0, 0). A^(k, w) shows a sharply defined quasi- 
particle band just below and at B = £?L2, see Fig. [5] (a) 



and (b). The fact that the residue Z^ does not vanish 
exactly at B = B\j± can be related to the finite tempera- 
ture. Also, we note that the single-particle residue is not 
fully converged in the intermediate field range, even at 
the lowest temperatures. 

In the static MF scenario, the two Lifshitz transitions are 
naturally present. As shown in Fig. [3j the single-particle 
residue Z^, calculated from the MF coherence factors at 
the Fermi energy, displays the expected step-like behav- 



The Lifshitz transition at B — £?L2 equally marks the 
transition from heavy to light fermions which is reflected 
in the steep increase of Z a as the magnetic field is ramped 
up further, see Fig. |6^b). This is in accordance with the 
notion of adiabatic continuity to free fermions which is 
expected in the limit of high magnetic fields, i.e. weak 
couplingPl Based on the fit = 200 DMFT results, we 
conclude that a continuous transition from low to high 
magnetic fields occurs, at least at and above this temper- 
ature. 



V. BEYOND DMFT 

The DCA calculates the k-dependent self-energy 
Y^ CA (uj n , K). This leads to the estimate for the residue 



?dca 



lim 

T->0 



(M(k f )) 



(15) 



0J n 



The map function M : kj — > K maps the Fermi momen- 
tum to the matching reciprocal patch. 
The 4-orbital DCA measurements agree with the 2- 
orbital DMFT results in the limits of strong coupling 
(small magnetic field) and weak coupling (large magnetic 
fields), see Fig.[7|[a). In the intermediate regime, around 
B = £>L2, deviations are detected in the magnetization 
as well as in the single-particle residues. The inclusion 
of spatial fluctuations softens the transition considerably. 
This can be understood from the notion of an effective 
Lande-factor gf which becomes lower when spatial fluctu- 
ations are present, since, on the two-site cluster, the local 
moment can be quenched not only dynamically but also 
via local singlet formation. The single-particle residue 
in the down-spin projection displays no sign of vanishing 
across the MMT (Fig. [^b)). 



VI. DISCUSSION 

Lifshitz transitions are continuous quantum phase 
transition which do not change symmetry but Fermi sur- 
face topology.^ Strictly speaking, they are defined for 
free fcrmion systems at zero temperature. Due to the 
unambiguous presence of quasiparticles, the notion of Lif- 
shitz transition can be carried over to the KLM. Driven 
by the external magnetic field, two consecutive Lifshitz 
transitions take place, at B = -Bli, 2, and the second one 
is identified with the MMT. This scenario is maintained 
when the /-moments are allowed to couple more strongly 
to the field by altering the ratio gf/g c - 
Collective effects challenging the quasiparticle coherence 
seem to be of minor importance during the MMT, even 
when B ~ r co h. Naturally, our calculation scheme 
is limited to the dominantly paramagnetic regime of 
the KLM. The choice of parameters, n c — 0.9 and 
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FIG. 7. (Color online) Magnetization M and single-particle 
residues Z CT from two-orbital (DMFT) calculations and four- 
orbital (DCA) calculations. The fc-vectors Ki = (0, 0) and 
K2 = (it, 7r) denote the relevant DCA patches. 



J/t > 1.3 place our res ults unambiguously in the para- 
magnetic phase pH 35 | 36 | pi rs t steps (Fig. n\ in a system- 
atic DCA study of larger clusters that can take into ac- 
count the Ruderman-Kittel-Kasuya-Yosida (RKKY) in- 
teraction between local moments leave the Lifshitz sce- 
nario at the MMT invariant. This is consistent with the 
fact that temporal fluctuations that generate the Kondo 
effect dominate the physics at the MMT. Close to a criti- 
cal point where the range of spatial fluctuations becomes 
large our approximation will fail and another modeling 
will be required. 

Transport signatures of the Lifshitz transition can be 
calculated with the Boltzmann transport approximation. 
Topological changes of the bands that cross the Fermi en- 
ergy can strongly influence transport measurements, in 
particular when these bands are shallow. This offers an 
explanation for the anomalies observed in Zeeman driven 
heavy-fermion systems!^ 

Compared to our results for the magnetic field dependent 
single-particle spectrum (Fig. [5]), similar results have 
been obtained for the ferromagnetic phase of the Kondo 
lattice model without external field termsP^ There, the 
spin-dependent shift of the quasiparticle weight is gen- 
erated dynamically and leads to the notion of a spin- 



selective Kondo insulating phase. 

Our results are applicable to heavy fermion compounds 
that have a magnetic field-driven Lifshitz transition at 
the coherence scale. 

The materials CeTiGeP and CeR^SiJ^ have a MMT 
at magnetic energy scales that are consistent with their 
estimated coherence temperatures. In our model, the 
critical metamagnetic field corresponds to the second Lif- 
shitz transition at £?L2- In this mechanism of competing 
energy scales we expect that the details of the band struc- 
ture are of secondary importance. This is in contrast to 
Lifshitz transitions at magnetic fields much below the co- 
herence scale where the details of the band structure are 
essential.^ 

The metamagnetic signatures of our model (Fig. [lj 
are similar to recent experimental data of the param- 
agnetic 4/-based compound CeTiGe which exhibits a 
pronounced first-order MMTP Its anticipated coherence 
scale, T co h ~ 55 K, is of the same order as the criti- 
cal magnetic field of /ijBmmt — 12.5 T, assuming in 
our model a g- factor gf « 7. Equally at lower fields, the 
magnetization is found to slightly change its slope, which 
might correspond to a first Lifshitz transition which in 
our model happens at -Bli- The experimentally observed 
distinct drop of the effective quasiparticle mass is in ac- 
cordance with our findings for the KLM (see Sec. IV I. 



Importantly, we find the MMT to be continuous both in 
the two-orbital DMFT and in the four-orbital DCA cal- 
culations and on the temperature scales we can access. 
CeRu2Si2 exhibits a continuous MMT and simultane- 
ously a Z eeman-driven topology change of the Fermi 
surface^. The magnetization increases seemingly linear 
as the magnetic field is increased towards the metamag- 
netic fields The critical field /j,qB mmt = 7.8 T matches 
the coherence temperature of T co h w 20 KP^2I when the 
g- factor in our model is assumed to be gf « 4. A Lifshitz 
transition at the coherence scale is therefore a plausible 
scenario for the MMT in CeRu 2 Si 2 . 



VII. CONCLUSION 

We have explored the Zeeman driven MMT in the 
Kondo lattice model which is considered to be the 
paradigmatic low energy model for heavy fermion sys- 
tems. Results for the paramagnetic metallic phase of 
the KLM are obtained in the framework of DMFT/DCA 
which can exactly account for the Kondo effect. 
Upon scaling the relevant energy scales with the lattice 
coherence scale the collapse of the magnetization data to 
a universal curve is observed, independent of the Kondo 
interaction. This data collapse has been confirmed for 
a range of Kondo couplings, temperatures and ratios of 
Lande factors. The pseudo spin nature of the /-orbitals, 
resulting from a Kramer's doublet, can be taken in ac- 
count with an effective Lande factor gf and the compe- 
tition of magnetic scale and coherence scale is invariant 
on the choice of gf. 
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We have traced the single-particle residue from low to 
high magnetic fields and report that it is continuous 
at the lowest temperatures our simulation can access. 
Two consecutive Lifshitz transitions occur as the field 
is ramped up and cause the change in topology of the 
spin-projected Fermi surfaces. This lead us to the find- 
ing that the MMT in the KLM is coincident with a con- 
tinuous Lifshitz transition. The absence of a singularity 
in the single-particle residue at the MMT excludes the 
Kondo breakdown scenario. 

At the temperature scale we can access, the sharp in- 
crease of magnetization at the MMT can well be ex- 
plained as a consequence of a continuous Lifshitz tran- 
sition in heavy fermion model systems where the Lande 
factor of the local spins is larger than the one for the 
itinerant electrons. In the course of this transition the 
excitations change their character from heavy fermions 
to light fermions. 

The recently observed first order nature^ of the meta- 



magnetic phase transition at a temperature T <C T co h 
remains an open issue. Of particular importance is un- 
derstanding if the KLM itself can account for the low 
temperature first order nature of the transition or if other 
competing energy scales such as coupling to the lattice^ 
have to be taken into consideration. 
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